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A symbolic method for /c-statistics^ 

E. Di Nardo, D. Senato* 



Abstract 

l/") , Trough the classical umbral calculus, we provide new, compact and easy to handle expressions of 

■ fc-statistics, and more in general of [/-statistics. In addition such a symbolic method can be naturally 

extended to multivariate case and to generalized fc-statistics. 

(N ■ 
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1 Introduction 

, In 1929, Fisher £|] introduced the fc-statistics as new symmetric functions of the random sample. The aim 

, of Fisher was to estimate the cumulants by free-distribution methods without using moment estimators. 

pH , He used only combinatorial techniques, fc-statistics are related to the power sum symmetric polynomials 

whose variables are the random variables (r.v.'s) of the sample. The Fisher point of view is described 
with a wealth of details by Kendall and Stuart in [Hj . The method is straightforward enough, however his 
execution leads to some intricate computations and some cumbersome expressions, except in very simple 
, cases. This is why many authors tried to simplify the matter later. 

CO . One of the most relevant contributions has been given by Speed [Jj in the Eighty. Speed resumed the 

, Doubilet approach to symmetric functions exploiting symmetric functions labelled by partition of a 

set rather than by partition of an integer. The entries in the transition matrices are computed via Mccbius 
, function generalizing and simplifying the presentation of fc-statistics theory. Nevertheless, in order to 



extend such theory to generalized fc— statistics, Speed must resort the tensor approach introduced by 
Kaplan in 1952. 

We follow a different point of view by using the high computational potential of the classical umbral 
calculus. This symbolic method was introduced by Rota and Taylor in 1994 (Hj and further developed in 
PP and [2]. From a combinatorial perspective, we revisit the Fisher theory as exposed by Kendall and 
Stuart, taking into account the Doubilet approach to symmetric functions. The umbral calculus offers a 
nimble syntax method that allows both the computation without using the Mcebius function and a natural 
extension to the multivariate case without bringing the tensor device into. What is more, this language 
clarifies the role played by the power sum symmetric polynomials in the expressions of fc-statistics. 

After recalling in Section 2 the strictly necessary umbral background, in Section 3 we put in umbral 
setting the four classical bases of symmetric polynomial algebra. Then we proceed in defining a general 
procedure to write down [/-statistics. Such a procedure is given in full details in Section 5, where analogous 
formulae for fc-statistics, /i-statistics and multivariate fc-statistics are given. 
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2 The umbral calculus language 

We start presenting the formal setting of the umbral calculus as introduced by Rota and Taylor in [Hj and 
further developed in [Q and [2]. We shall confine our exposition to what is necessary to the aims of this 
paper. 

The umbral calculus is a syntax consisting of the following data: an alphabet A = {a, (3, . . .} whose 
elements are named umbrae; a commutative integral domain R whose quotient field is of characteristic 
zero; a linear functional E, called evaluation, defined on the polynomial ring R[A] and taking values in 
R and such that E[l] = 1, E[a { f3 j • • • 7 fe ] = E[c^]E[p j \ ■ ■ ■ E[y k ] for any set of distinct umbr ae in A and 
for i,j,...,k non-negative integers (uncorrelation property); an element e G A, called augmentation, such 
that E[e n ) = for n > 1; an element u £ A, called unity umbra, such that = 1, for n > 1. 

The support of an umbral polynomial p G R[A] is the set of all umbrae occurring in p. Two umbral 
polynomials are said to be uncorrelated when their supports are disjoint. 

An umbra carries the structure of a r.v., while making no reference to a probability space if the 
evaluation E is considered as the expectation operator. The moments of an umbra a are the elements 
a n G R such that £?[a n ] = a n for n > 1. We said that the umbra a represents the sequence 1, oi, 0,2, .... 
The singleton umbra x is the umbra whose moments are all zero, but the first equal to 1. The factorial 
moments of an umbra a are the elements ar n -\ G R corresponding to the umbral polynomials (a) n = 
a(a — 1) • • • (a — n + 1), n > 1 via the evaluation E, i.e. £?[(a)„] = a,( n ). The Bell umbra j3 is the umbra 
such that E[(P) n ] = 1, n > 1. The umbra (3 and the umbra x allow respectively a symbolic tool in order to 
handling composition and inversion of formal power series (for a detailed exposition of such two umbrae 
with their properties see £Q and [2]). 

Two umbrae a and 7 are said to be similar when = -£[7™] for all n > 1, and we will write 

a = 7. Then two similar umbrae represent the same moment sequence. Given a sequence 1, a\, 02, ... in 
R this is represented by infinitely many distinct and thus similar umbrae. Such a circumstance permits 
to deal with moment products by means of the uncorrelation property. Indeed, it is a^a, 7^ i?[a I a J ] with 
a- L = E[a l ] and aj = E[a J ], as well as a;Oj = E[a l a 3 ] with a = a and a uncorrelated with a. So, given 
n G N, the sequence 



Let p and q be two umbral polynomials. We said that p is umbrally equivalent to q iff E[p] = E[q], in 
symbols p ~ q. This last equivalence relation turns out to be very useful in defining and handling umbra 
generating functions. We point out that all operations among umbrae correspond to analogous operation 
in the algebra of generating functions. Nevertheless, in the following we make no mention of generating 
functions (referring jHj to a formal exposition). 

Thanks to the introduction of similarity notion, it is possible to extend the alphabet A with the so- 
called auxiliary umbrae derived from operations among similar umbrae. This leads to the construction of 
a saturated umbral calculus in which the auxiliary umbrae are treated as elements of the alphabet (cf. [H])- 
Let a , a , . . . , a be n uncorrelated umbrae similar to an umbra a. The symbol n.a denotes an auxiliary 
umbra similar to the sum a + a + • • • + a . The symbol or" denotes an auxiliary umbra similar to the 
product a a ■ ■ ■ a . Properties of such auxiliary umbrae are extensively described in pQ and they will 




gives the moments of a + a . 



be recalled whenever it is necessary. We will assume disjoint both the support of n.a,m.a and a n ,cc m 
whenever n m. If p and q are correlated umbral polynomials, then n.p ~ pi + • • • + p n is correlated 
to n.q ~ qi + ■ ■ ■ + q n , and pt is correlated to qi but uncorrelated to qj with i ^ j. In pp, the following 
identity is stated: 

i 

E[(n.a) 1 ] = ^2(n) k B i>k (ai,a2, . . . ,a.j_ fc+ i) i > 1, (1) 
fc=i 

where -B^fc are the (incomplete) exponential Bell polynomials and otj is the i-th moment of a. Moreover , 
it is easy to verify that E[(a- n ) 1 ] = a" for i > 0. 

Two umbrae a and 7 are said to be inverse to each other when a + 7 = e. The inverse of the umbra 
a is denoted by — l.a. Note that, in dealing with a saturated umbral calculus, the inverse of an umbra is 
not unique, but any two inverse umbrae of the same umbra are similar. 

Replacing the integer n in n.a with an umbra 7, we obtain the auxiliary umbra 7. a whose moments 

are 

i 

E [(j- a Y] = ^2g(k)B t ,k(ai,a2, ■ ■ ■ ,Oi-k+i) i > 1, (2) 
fc=i 

where are the factorial moments of 7. In particular f3.a is called a— partition umbra and its moments 
are the (complete) exponential Bell polynomials (cf. [Q). Moreover x- a is called a— cumulant umbra and 
a.x is called a— factorial umbra, with moments equal to the factorial moments of a (cf. [2])- In particular 
it is 

p. X = u = X .f3. (3) 

Again, replacing the umbra 7 in j.a with the umbra 7./?, we obtain the composition umbra of a and 7, 
i.e. 7. /?. a. The compositional inverse of an umbra a is the umbra a < ~ 1> such that a <_1> ./3.a = x = 
a./?.a <_1> . In particular it is 

w < - 1> ./3 = X = m <_1> -/?, (4) 

where m <_1> denotes the compositional inverse of u. Via the umbral Lagrange inversion formula (cf. U|), 
the moments of u < - 1> are E[(u < - 1> ) n ] = (-l)"" 1 ^ - 1)!. Finally it is 

X.X = « <_1> - X-(-x) = (-«) < - 1> - (5) 

The disjoint sum of a and 7 is the umbra whose moments are the sum of n-th moments a n and g„ of a 
and 7 respectively, in symbols (a+7)" ~ a" + 7™ (cf. [2]). For instance it turns out 

X-u+X-1 = X-(<* + i) ( 6 ) 

that is the well known additive property of cumulants. In the following, we denote by +„a the disjoint 
sum of n times the umbra a. 



3 Umbral symmetric polynomials 

A partition of an integer m is a sequence A = (Ai, A2, . . . , Xt ), where Xi are weakly decreasing and 
J2t=i = m - The integers Ai are said parts of A. A different notation is A = (l ri , 2 r2 , . . .), where r 2 ; is the 
number of parts of A equal to i. The monomial symmetric polynomials in the variables ct\, 012, ■ ■ ■ ,a n are 
m\ = c^i 1 ' ' ' c^t 1 ■> where the sum is over all distinct monomials having exponents Ai, . . . , At. When A 
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ranges in the set of partition of the integer to, m\ is a bases for the algebra of the symmetric polynomi- 
als. There are different other bases; we recall just those necessary in the following: the r-th power sum 
symmetric polynomials s r = Yl^—j the fc-th elementary symmetric polynomials = 
where the sum is over 1 < j\ < j% < ■ ■ ■ < jk < n, and the r-th complete homogeneous symmetric 
polynomials hi = X)|A|=i mx - 

When the umbrae £*!,••• , a n are uncorrelated and similar to each other, these four classical bases of 
the symmetric polynomial algebra can be represented by means of the umbral polynomial n.(xot l ) and 
its moments. Propositions 1, 3 and 4 are stated under this hypothesis. For the power sum symmetric 
polynomials it is 

s r ~ n.a r ~ n.(xa r ). 
Note that s r are umbrally equivalent to the moments of + n « (cf. 0)- 
Proposition 1 (Umbral elementary polynomials) 

Proof. For k = l,...,n the result follows applying the evaluation E to the multinomial expansion of 
[n.(xa)] fe ~ (xi a 'i + ' ' '+Xn a n) k and observing that terms having powers of x greater than 1 vanish. So just 
A:! monomials of the form Xji a ji Xh a h ' ' ' Xjk a jk nas an evaluation not zero. Instead for k = n+1, n+2, . . . 
the result follows observing that at least one power of x greater than 1 occurs in each monomial of the 
multinomial expansion. ■ 

Proposition 2 It is 

X-n.(xa) = u <_1> (+„£*) n.( X a) = p. [u^^ (+„«)] (8) 

Proof. The first equivalence in JBJ) follows from © observing that 

X-n.(xa) = X-ixmi + ■ + X n a n ) = +nX-(x«) = +„M <_1> a. 

The second one follows from the first, making the right dot-product with (3 and recalling J2J). ■ 

Equations JHJ) are the umbral version of the well-known relations between power sum symmetric polyno- 
mials s r and elementary symmetric polynomials e^. Indeed s r is umbrally equivalent to the moments of 
(+ n a) and is umbrally equivalent to the moments of n.(x a )- The umbral expression of m\ requires the 
introduction of augmented monomial symmetric polynomials fri\. Let A = (l n , 2 r2 , . . .) be a partition of 
the integer to, such polynomials arc defined by 

™* = S a * " ■ a >i a ii +1 ■ ■ ■ a l 1+ . 2 ■ ■ ■ • ( 9 ) 

1 < h ■ ■ ^irj +i • • •jV-i+ra ■ ■ -<n 
The proof of the next proposition follows the same patterns of Proposition ^ 

Proposition 3 It is 

™A ~ [n.(xa)] ri [n.(xa 2 )} T2 
The next corollary follows recalling that ni\ = toa/^i!^! • • •]. 



5 



Corollary 3.1 (Umbral monomial polynomials) 

[n.( X a)r Mx<* 2 )] ra 



to a 



Proposition 4 (Umbral complete polynomials) 

[-n.(- X a)} m ~ m\h m to =1,2,.... (10) 

Proof. It is 



-ii 



(- X a)] m ~ [_i.(- Xiai ) + ... + -l.(-x„a n )]' 



m! ^ [-l.(- X ')]- n ([-l-(-x")] 2 )- r2 ~ 



|A|=m V ' V ; 

and the result follows from Corollary 13. II being ([— 1.(— x)D' n — 

Note that equivalences (0) and (|1C)|) are the umbral version of the well-known identities: 



n 1 

= IK 1 E**** ^ ( } . 

fe i=l fe lll=lV 1 / 



Proposition 5 /£ is 

-X.n.(-xa) = (-«) <_1> M - = P- [(-«) <-1> (+„«)] (H) 

Proof. The results follow from ijHJ) replacing the umbra \ with — \ and recalling that u <_1> = x-X m ust 
be replaced with (— it) <_1> = — x( — X)- ■ 

Equations (|ll|l are the umbral version of the well-known relations between power sum symmetric poly- 
nomials s r and complete symmetric polynomials hk- Indeed s r is umbrally equivalent to the moments of 
(+ n oi) and hk is umbrally equivalent to the moments of — n.(— X a )- 

4 [/-statistics 

Let Xl, X2, ■ ■ ■ , X„ be n independent r.v.'s. A statistic of the form 

u = -^- k ^( x h,x j2 ,...,x jk ) 

where the sum ranges in the set of all permutations (ji,j2, • • • ,jk) of k integers, 1 < ji < n, is called 
[/—statistic |5J- If Xi, X%, ■ ■ ■ , X n have the same cumulative distribution function F(x), U is an unbiased 
estimator of the population character 0(F) = f ■■■J $(21, ■ ■ • ,Xk)dF(xi) ■ ■ -dF(xk). In this case, the 
function $ may be assumed to be a symmetric function of its arguments. Often, in the applications, $ 
is a polynomial in Xj's so that the [/—statistic is a symmetric polynomial. By virtue of the fundamental 
theorem on symmetric polynomials, the [/—statistic can be considered a polynomial in the elementary 
symmetric polynomials. The following theorem is an umbral reformulation of the above statement. 

Theorem 4.1 ([/—statistic) If A = (l ri , 2 1 " 2 , . . . , ) is a partition of the integer to < n then 

K-J^r • • • ^ ? ^ r [n.(xa)PMxa 2 )P ■ ■ ■ (12) 
where ji £ {1, 2, . . . , n} and Y] - rj = k. 



G 

Proof. The result follows observing that [n.{xot l )Y i — {n-x) ri {a l )' ri and (n.\) k ~ (n)k- ■ 

The formula (|12|) states how to estimate moment products by using only n information drawn out the 
population. Then the symmetric polynomial on the right side of Ijl2|) is the U— statistic of the uncorrelated 
and similar umbrae on, o>2, . . ■ ,a n associated to (a^)^ 1 ( Q j 2 )' r2 ' • • ■ 

Example 4.1 Moment powers. Set r\ = 2 and k = 2, from (|12|l the U— statistic associated to or 2 ~ a\ is 

a 2 ~ — ; r\n.{Ya)} 2 ~ — : Vttja,-, n > 2. 

n{n-\y yA n n(n-l)^ J 

Example 4.2 h-statistics. As it has been shown in [2|, it is 

where a Ql is the central umbra of a about a\ = E[a] and a = a is an umbra uncorrelated with a. 
Replacing the product cc k (a ) r ~ k with the corresponding U— statistic I) 1211 . it results 

KT ~ £ ( £ ) (ziH[ n . (xa )]V (xa ^). (13) 
fc=o ^ ' \ n ) k+1 

When the products [n.(xc t )] k n.(xct r ~ k ) are expressed in terms of power sum symmetric polynomials, we 
have the so called /i-statistics. On this matter, we will give more details in the next section. 

5 ^-statistics 

The n-th fc-statistic k n is the unique symmetric unbiased estimator of the n-th cumulant K n of a given 
statistical distribution, i.e. E[k n ] = n n . The /c-statistics can be expressed in terms of the sums of the r-th 
powers of the data points. In this section we give an umbral syntax that provides a general computational 
method to generate such expressions. To this aim, we digress to introduce the exponential Bell umbral 
polynomials. 

The most widespread expression of incomplete exponential Bell polynomials is referred to partition 
of an integer. Of course, it is also possible to express such polynomials referring them to partition 
of a set. Here we follow this last point of view. We denote by Hi.k the set of all partitions of the set 
[i] = {1, 2, . . . ,i} in k blocks. Let tt = {Ai, A2, . . . , A^} be an element of n^.. Then it is Bi^{a\, 02, . . .) = 
Svren fc a ni a„ 2 • ■ ■ a nk where \Aj \ = rij,j = 1, 2, . . . , k, as we will suppose from now on. Let us consider 
the following symmetric umbral polynomial: 

Hj.jfcCax.aa,--- ,a k )= ]T <^-aJ», (14) 

where ai, ot2, ■ ■ ■ ,o>k are uncorrelated umbrae similar to a and a nj = E[a n '],j = 1,2, ... ,k. Obviously it is 
E[Bi t k] — Bi t k, so that any expression containing the polynomials Bi ^ could be replaced with an umbrally 
equivalent expression containing the polynomials Bi t k- The polynomials Bi t k will be called (incomplete) 
umbral exponential Bell polynomials. The combinatorics underlain the polynomial Bi t k is the following: 
the set [i] is partitioned in k blocks, to each of them one associates the umbra a nj obtained firstly replacing 
the elements in the j-th block with the umbra a and then labelling all blocks so that powers belonging 
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to different blocks result uncorrelated. Replacing in (fbfjl the products a™ 1 a^ 2 • • • a^ k with the umbrally 
corresponding U— statistic lf[2")) . we get for i < n 

B hk {a 1 ,a 2 ,--- ,a k )~ ^ -^—n.{ X a ni )n.(xa n2 )---n.(xa nk ), (15) 
nen,, k [n)k 

by which we are able to give the umbral fc-statistics. Indeed, the a~ cumulant umbra n a is similar to X- a i 
so that, being (x)k — (u < ~ 1> ) k , from J5J) and l(T5)l it is 

fra)*-^-!)*-*!*^! ^ n.( X a">.(x^) •••n.( X a"*). (16) 

Since n.(xa ni ) is umbrally equivalent to a symmetric power sum polynomial, equation (|16f) gives the 
moments of the a— cumulant umbra in terms of power sum polynomials, i.e. the umbral form of k- 
statistics. Note that the symmetric power sum polynomials in l|16|) are correlated. So, in order to make 
formula (|16f) effective, we need a device by which to evaluate the product n.^a" 1 ) n.(xa™ 2 ) • • • n.(xct 7lk ). 
To this aim, by using the umbral exponential Bell polynomials (11411 . the moments of the umbra n.(x a ) 
can be evaluate from J2J) recalling n.a ni ~ (+ n a) ni and the second equivalence in (JSJ 

i 

[n.{ X a)} l ~Yl 12 (u' < - 1> ) n Hw'' < - 1> r 2 ---(/' < - 1> ) nfc [n.a mi ][n.a' m2 ]---[n.a"'" fe ]. (17) 
fe=i u-en iifc 

The previous umbral equivalence is suitable to be generalized to the product of umbral polynomials 
[n.(xpi)] [n.(x?>2)] • • • [i-(xPi)] with no disjoint support. Indeed, it results 

i 

[n.( X Pi)]---[n.(m)]^J2 E (u < - 1> r---(u"<- 1 >r k [n.p' Ai ] ■ ■ ■ [n.P' Ak ] (18) 

k=l Trelli, k 

where P Aj = Yit= l Pjt an d Pjt are the polynomials indexed by the elements of the block Aj, as we will 
suppose from now on. Equivalence l(T%|) is the device required to generate fc-statistics. Indeed, setting 
i = k and p t = a nt for t = 1,2, ■ • ■ ,i, it results: 

k 

[n.{xa ni )}---[n.( X a nk )} ^ {u^-^)" 1 ■ ■ ■ (u" < - 1> ) n * [n.a'" 11 ] ■ ■ ■ [n.a"™'] (19) 

where = 52™= l n jt an d n jt are indexed by the elements of the block Aj. Note that the power sum 
polynomials on the right side of l|19|) are now uncorrelated so such equivalence gives augmented monomial 
symmetric polynomials in terms of power sum polynomials, translating Kendall and Stuart tables read 
downwards |Hj. Instead the following formula leads to a symbolic translation of the Kendall and Stuart 
tables read across (i.e. power sum polynomials in terms of augmented monomial ones): 

i i 

(n.pi)...(n.pi)^£(n) fc £ P Al ■ ■ ■ P Ak ~ ]T ]T n.( X P Al ) • ■ ■ n.{ X P Ak ). (20) 

fe=i Tren^fc fc=i7rGn iifc 

The first equivalence in (|20J) is obtained from 0) trough analogous considerations used to state (|18|l : the 
second equivalence comes from llll't replacing o^ 1 with the umbral polynomial pi, i.e. 

PM P A 2 --- P A k - T^n.( X PA 1 )n.(xPA 2 )---n.( X P Ak ). 



Set in Pt = OL nt for t = 1, 2, i, we have 

i 

(n.a ni )- ..(n.a n< ) Z (™-Xa mi )("-Xa m2 ) • • • (™-Xa mfc ) (21) 

fc=i 7rgn iiti 

where TOj = Y^t=i n jt an( i n jt arc indexed by the elements of the block Aj. 

Example 5.1 h-statistics. In l|19|l set ri\ = n — k and n 2 = • • • = Jifc+i = 1; from (|13|) we get the umbral 
expression of /i-statistics. 

Example 5.2 Joint cumulants. Let Pi,p%, ■ ■ ■ ,Pi be umbral polynomials. Replacing n with \ i n the first 
equivalence in l|2()[l . we have 

i 

(X-Pi)(X-P2)---(x-Pi)-J2W* E P M--- p A k - (22) 
fc=i 7ren ifc 

When the umbral polynomials pi are interpreted as r.v.'s, equivalence 11221) gives their joint cumulants. 
So we will call (x-Pi)(X-P2) • ■ ■ (x-Pi) the joint cumulant of pi, . . . ,pi. Note that, setting p t = a for t = 

1, 2, . . . , i, one has the i-th ordinary cumulant (x- a Y ■ Through this equivalence it results X-{Pi H ^Pi) = 

X-Pi + ' ■ '+X-Pi - Now suppose to split the set {pi,p 2 , ■ ■ ■ ,Pi} in two subsets {pj 1 , . . . ,pj t } and {pk x , ■ ■ ■ ,Pk s } 
with s + 1 = i, such that polynomials belonging to different subsets are uncorrclated. Then it is 

(X-Pi)(x-P2)---(X-Pi)^0. (23) 

Indeed, setting P = J2i=iPji an( i Q = J2i=iPkn sucn polynomials are uncorrelated, so that x-(P + Q) = 
X-P+X-Q from ©. Equivalence lj2"3"|l follows observing that, due to the disjoint sum, products involving 
powers of x-P and x-Q vanish. When the umbral polynomials pi are interpreted as r.v.'s, equivalence ()23j) 
states the following well-known result: if some of the r.v.'s are uncorrclated of all others, then their joint 
cumulant is zero. 

Example 5.3 Multivariate k-statistics. Equivalence 1)20(1 allows a compact expression of multivariate k- 
statistics. Replacing n with x in its second equivalence, we construct the U— statistic of the joint cumulant 

(X-ft)(X-fc)-(X-ft)-trr E n.( X P Al )n.(xP A2 )---n.(xPA k ). (24) 

Again, in the product on the right side of l|24|l the umbral polynomials are correlated. In order to make 
effective the computation, it is necessary to rewrite (|24|) by using equivalence (|18[) with Pa ± instead of pt- 
For instance, in order to express fc 2 i, set in <|24|) i = 3 and p\ = P2 — cxi,ps = 012- It results 

(X-"i) 2 (X-a2) ^ -n.(xala 2 ) + {2n.(xa 1 )n.(xaia 2 ) + n.(xc4)n.(xa 2 )} 



{[n.( X ai)] 2 n.( X a 2 )} (25) 



Set s P: q ~ n.(af a 2 ). It is 



n.(x«i) n.(xaia 2 ) : 


, („<- 


" 1> ) 2 n. 


n.(xa 2 )n.(xa 2 ) = 


- (u<" 


- 1> ) 2 n. 


{n.(xa 1 )} 2 n.(xa 2 ) 


- (u <_ 


- 1> ) 3 n. 




1- u 


_1> (u< 



) 2 n.(a 2 a 2 ) + (w <_1> )' 2 n.a' 1 n.(aia: 2 ) ~ — s 2 ,i + si,o Si,i (26) 
) 2 n.(a 2 a 2 ) + (w <_1> )- 2 n,.a' 2 n.a 2 ~ — s 2 ,i + s 2l o So,i (27) 
) 3 7i.(a 2 a 2 ) + (u < ~ > ) n.a'i n.a.\ n.a 2 
) 2 [n.a 2 n.a 2 + 2n.a 1 n.(aia 2 )] 
2 S2,i - s 2 ,o s ,i - 2 81,0 si,i + s 2 s o,i (28) 
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Equivalence H'26J) comes from (|18|l setting i = 2, pi = cti andp 2 = 0*10:2, equivalence (|27|l comes from ()18(l 
setting i = 2, p\ = a\ and P2 = o? 2 , equivalence l|28|) comes from (11 81) setting i = 3, p\ = P2 = Qi and 
P3 = 012- Substituting the above equivalences in 125fl and rearranging the terms, we have the expression 
offc 2 i, 

fc 2 i ~ (x.ai) 2 (x.a 2 ) ^ 7-tt- [n 2 s 2 ,i - 2nsi. si.i - ns 2 . s ,i + 2 s 2 s ,i] . 
W3 

The expression of generalized fc— statistics (as well as the multivariate ones) in terms of power sums comes 
from (|22J) replacing some of the umbrae x with uncorrclated ones and then constructing the corresponding 
[/—statistics. 
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